Submitted to Ap. J. 



in 
o 
o 

(N 
u 

< 



OBSERVATIONAL EVIDENCE OF THE FORMATION OF 
CYANOPOLYYNES IN CRL618 THROUGH THE 
POLIMERIZATION OF HCN 

Juan R. Pardo, Jose Cernicharo & Javier R. Goicoechea 1 

IEM - Departamento de Astrojisica Molecular e Infrarroja, CSIC, Serrano 121, E-28006 

Madrid, Spain. 

1 Present address: LERMA-Ecole Normale Superieure (UMR 8112 du CNRS), 24 Rue 

Lhomond, 75231 Paris, France. 



pardo , cerni , javier@damir . iem. csic . es 

in 

Csj ■ ABSTRACT 

o 
in 
o 



The abundance ratio of consecutive members of the cyanopolyynes family has 
been explored in CRL618 using data acquired in a complete line survey cover- 



ing the frequency range 81-356 GHz. The J up range explored for the different 
molecules is the following: 1 to 4 for HCN and HNC, 9 to 39 for HC 3 N, 31 to 
133 for HC 5 N, and 72 to 85 for HC 7 N (not detected beyond J up =85). The lowest 
vibrationally excited state of HC7N (z/ 15 at 62 cm -1 ) has been tentatively de- 
tected. Data analysis has been performed by extending our previous geometrical 
and radiative transfer model of the slowly expanding envelope (SEE) surrounding 
the compact central continuum source of CRL 618, that was established from the 
study of rotational lines in several vibrationally excited states of HC3N. The new 
lines analyzed here require to model the high velocity wind (HVW) component 
and the colder circumstellar gas, remnant of the AGB phase of CRL618. The de- 
rived HC3N/HC5N and HC5N/HC7N abundance ratios from this set of uniformly 
calibrated lines are between ~ 3 and 6 in the different regions, similar to standard 
values in the CSM and ISM, and consistent with previous estimates obtained from 
ISO observations and chemical models. However, the abundance ratios of HC3N, 
HC5N and HC7N with respect to HCN are at least two orders of magnitude larger 
than those typical for AGB C-rich stars, such as IRC+10216. This fact indicates 
that, in the short transition toward the Planetary Nebula phase, HCN is quickly 
reprocessed into longer cyanopolyyne chains. A similar behavior was previously 
found in this object for the polyacetylenic chains (C 2n H 2 ). 
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Subject headings: stars: post-AGB-stars: carbon-rich - stars: circumstellar mat- 
ter - stars: individual: CRL618 - ISM: molecules - radio lines: stars 



1. Introduction 

CRL 618 is a C-rich protoplanetary nebula (PPNe). Its chemical richness and its com- 
plex morphology with a central BO star surrounded by an ultracompact HII region, and a 
thick molecular and dusty envelope (Bujarrabal et al. 1988) with different outflows (Cer- 
nicharo et al. 1989, hereafter CER89), are now well known thanks to many detailed obser- 
vational studies. Detailed chemical models for this object have been developed (Cernicharo 
2004, hereafter CER04) indicating the rapid evolution of the central star and its influence 
on the circumstellar ejected material. Changes seem to occur at short -almost human- time 
scales. The most recent observational results on the morphology are the interferometric 
observations presented in Sanchez-Contreras & Sahai, 2004. Concerning the chemical com- 
position, a line survey from 80 to 275 GHz has been recently completed (Cernicharo et al. 
2005). 

After the study of the physical conditions in the slowly expanding envelope (SEE) around 
the ultracompact HII region of CRL 618 using vibrationally excited states of HC 3 N (Pardo 
et al. 2004, hereafter PAR04), we focus now on the entire cyanopolyynes family: HCN, 
HNC, HC 3 N, HC 5 N and HC 7 N (longer members have not been detected in our survey). 

HC 5 N is the second molecule, after HC 3 N, in terms of total number of lines detected 
in the millimeter spectrum of CRL 618. It is a linear molecule with a rotational constant 
of 1331.330 MHz and a dipole moment of 4.33 Debyes, first discovered in space by Little et 
al. (1978). The two lowest energy doubly degenerated bending modes (vu and u 10 ) have an 
energy of 105 and 230 cm" 1 respectively (Yamada et al., 2004). On the other hand, HC7N was 
first discovered in space by Kroto et al. (1978). Its rotational constant and dipole moment 
are 564.0011 MHz and 4.82 Debyes respectively. HC3N/HC5N and HC 5 N/HC 7 N abundance 
ratios between 1.4 to 3.0 were found in molecular clouds since the earliest studies (Snell 
et al. 1981). Finally, HCN, with a dipole moment and rotational constant of 2.98 Debyes 
and 44315.976 MHz respectively, is one of the most abundant molecules in the interstellar 
and circumstellar media. In many cases HCN lines are optically thick so that observations 
of isotopologues or vibrationally excited states are preferred to probe molecular clouds and 
circumstellar envelopes. The other two molecules discussed in this paper have the 
following dipole moments and rotational constants: 3.73 Debyes and 4549.059 
MHz (HC 3 N), and 2.98 Debyes and 45331.7845 MHz (HNC). 
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The goal of the study conducted in this paper is to get abundance ratios of the different 
members of the cyanopolyynes family in CRL 618 and to compare them with those observed 
in prototypical AGB objects such as IRC+10216. In addition, by using new polar species 
we aim at improving and extending our previous model of CRL618 (see PAR04). This is 
another step toward a final model for the whole millimeter wave spectrum of this source. The 
observational procedure and SEE model developments were described in detail in PAR04, 
we only give here a short summary in sections 2 and 3. The analysis of the HCN, HNC, 
v—0 HC 3 N, HC 5 N, and HC 7 N data, with the extension of the model, is presented in section 
4. Comparisons of the results with chemical model predictions, and other discussions, are 
presented in section 5. The work is summarized in section 6. 



2. Observations 

The observations presented in this paper are part of two line surveys of CRL618. The 
first one, now complete, has been carried out with the IRAM-30m telescope from 1994 to 
2002 (Cernicharo et al. 2005, hereafter CER05) with a frequency coverage of 81-279 GHz 
except for frequencies with high atmospheric opacity around 119 and 183 GHz. The second 
one, still on-going but complete for the species we focus on in this work, is performed with 
the CSO telescope between 280 and 356 GHz. For both instruments the frequency 
resolution used has been around 1 MHz with ~0.5 GHz wide spectrometers. The 
J up range covered for rotational transitions of HCN and HNC is 1 to 4, for HC3N 9 to 39, 
for HC 5 N 31 to 133, and for HC7N it is 71 to 316, although the lines are not detected beyond 
J«p=85. See PAR04 for details on the observational procedures at the two facilities. 



3. Previous Model Summary 

A model has been developed to study the physical conditions that explain the whole ob- 
served millimeter wave spectrum of CRL618. The precise knowledge of the spectral behavior 
of the continuum emission (see PAR04) is paramount since it plays a very important role 
in the observed line profiles (specially in the emission-to- absorption line ratio). The model 
consists of a gas envelope around a central continuum source, considered spherical with size 
and effective temperature adjusted to fit the continuum IRAM-30m data (PAR04) and as- 
sumed to keep the same spectral index in the range 280-356 GHz (CSO observations). The 
description of the gas component allows different shells with different physical conditions to 
be included in the model. The geometry is not restricted to be spherical. Elongated shells 
with both radial and azimuthal velocity components can be defined. 
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Table 1: Observational parameters (velocity centroid, width and integrated area) obtained 
from gaussian fits of the absorption and emission components of selected HC5N lines in 
CRL618. See beginning of section 4 for more details. 
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Since the observations reveal that different gas regions are traced by different molecular 
species, we proceed in our analysis in different steps. The first one has been to study 
those species arising from the inner ~1.5" SEE, with outflow velocities of ~ 10 to 15 km-s -1 , 
turbulence velocities around 3.5 km-s" 1 , and temperatures in the range 250-275 K. The 
full description of the SEE model is given in PAR04. The extension of the 
SEE model to include both the colder and outer circumstellar gas (necessary to 
explain the rotational lines in the ground vibrational state of HC 3 N, HC 5 N and 
HC 7 N) and the high velocity wind component (seen in lines of some abundant 
species, including v =0 HC 3 N) is one of the targets of this work. This is discussed 
in detail in the next section. 



4. Analysis and Discussion 

The first step of the analysis has been to perform gaussian fits of selected HC 5 N lines 
in the 3, 2, and 1.3 mm windows. A summary for a set of lines representative of the 
different vibrational states and frequency ranges is presented in Table 1 where the 
areas are expressed in K-km-s -1 (T*a' v lsr), parameters that could not be determined are 
marked with a dash, and parameters that had to be fixed in order to obtain a reasonable 
fit are highlighted in bold face. Also, when the velocity, width and area of one component 
are left blank in the table it means that the component is not detected in the particular 
transition. These observational results already reveal an interesting behavior that is crucial 
for the analysis: lines in the v =0 state have absorption and emission components that are 
in general 15 to 25 km-s -1 wide at 3 and 2 mm, with the absorption centered at around 
-37 to -42 km-s -1 , whereas lines in vibrationally excited states, in the same spectral region, 
appear narrower (~ 7 to 15 km-s -1 ) with the absorption component (sometimes the only one 
detected) centered at ~ -28 to -33 km-s -1 . These differences tend to disappear as frequency 
increases and the line profiles in the v=0 state are dominated by the emission component 
only. These facts suggest that the rotational lines in the ground state of HC5N are mostly 
formed in a region with quite different physical conditions than the one responsible for the 
lines in vibrationally excited states. This is carefully explored in this section. 



4.1. Vibrationally excited HC 5 N 

The latest laboratory data on the v\\ and v\§ vibrational states of HC5N, their overtones 
and their combinations, have been presented by Yamada et al. (2004). In our line survey, 
we have detected rotational lines in the following vibrationally excited states of HC 5 N (most 
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of them for the first time in space): uu, 2un, v w , and Sun, at energies of 105, 210, 230, and 
315 cm -1 respectively. Since these lines also should arise from the SEE analyzed in PAR04, 
we have just tried to find the HC3N/HC5N abundance ratio that best fits the observations, 
keeping the same SEE physical parameters of PAR04. The result is HC3N/HC5N ~ 3±1 
(see Fig. 1). This result confirms previous estimates made from vibrational bands of both 
species seen in absorption by ISO around 15 /im (Cernicharo et al. 2001). Obviously, the fit 
cannot be perfect for so many lines. In particular, Fig. 1 shows that some absorption exists 
between -50 and -40 km-s -1 for J up =35 in the v u state that is not reproduced by the model. 
This can be due to a fraction of a higher velocity outflow that intersects the line of sight 
toward the continuum source that has not been considered in the model. This absorption 
component is otherwise clearly visible in the rotational lines inside the z/7 state of HC3N, 
but does not affect the other lines shown in Fig. 1. The detections in the HC5N vibrational 
states extend from J up =31 to J up =104 for i/ n , J up ~87 for 2z/ n , J up ~77 for v w , and J up ~53 
for 3z/n. The fact that for HC3N we detect rotational lines from vibrational states a factor 
3-4 higher in energy than for HC5N is consistent, for the sensitivity of the survey, with the 
derived HC3N/HC5N abundance ratio, and the dipole moment and partition function ratios 
of both species. 



4.2. 12 C/ 13 C ratio in the SEE 

The 12 C/ 13 C isotopic ratio in the SEE can be derived using HCC 13 CN, HC 13 CCN and 
H 13 CCCN in the u 7 state. The v$ state is also detected but the signal-to-noise ratio is too low 
for any estimation. Figure 2 shows the results of the SEE model for ( 12 C/ 13 C) = 15 (PAR04). 
We estimate an error bar for this value of ~ ±2. Kahane et al. (1992) gave lower limits 
for this ratio of 18 and 3.2 depending on which transition (J=l-0 or J=2-l) of both 12 CO 
and 13 CO was used. Our estimate gives a much better constraint, and restricted to the SEE 
region, because it relies upon tens of lines without the saturation problems of CO. 



4.3. Ground state HC 5 N 

Using the result found in 4.1, we have run the model for the detected rotational lines 
of HC 5 N in its ground vibrational state. From these calculations it becomes clear that 
the SEE alone cannot account for these lines. A significant contribution from gas at lower 
temperatures must exist. In the individual lines, there is no significant evidence of 
the high velocity wind (HVW) seen in the most abundant molecules (see below). Therefore, 
we have expanded the model presented in PAR04 to include a cold circumstellar shell (CCS), 
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created during the AGB phase (see figure 3). We have not simply considered a spherical 
geometry for this cold component in order to account for the optical/infrared images of the 
CSE that show that the envelope has a bipolar structure. This was also pointed out for the 
molecular gas in CER89 and Neri et al. (1992). The parameters that define the CCS in our 
model are: 

• Rccsi an d Rccs2'- Inner and outer boundaries of the shell (angles in arcseconds). 

• ar- Truncation angle (degrees), used to take into account the fact that this shell needs 
to be divided into at least two zones because the molecular density is probably larger 
in the polar regions along the axis of the high velocity wind. This fact is revealed by 
the ratio between emission and absorption in the profiles of the analyzed v=0 HC5N 
rotational lines. 

• ^ot S x'- Rotational temperature (K) for species X in the CCS. 

• [HC3N] (at R=Rccsi)'- Number density of HC3N at the inner part of the shell for 
a > ax (for a < ay the density will by multiplied by a factor p ccs ). HC 3 N will always 
be our reference molecule for the models of CRL 618. In each region we will consider 
a constant ratio with respect to it for each molecule. 

• decs'- The radial distribution of HC 3 N abundance is considered to vary as [HC 3 N^] = 
[RC 3 N Rccsi }(R/R CC si) dccs in the CCS. 

• v c*C5 an d v tur b,ccs'- radial and turbulence velocity of the gas in the CCS. 

Just as we did in PAR04 for the physical parameters of the SEE, we had to run a 
considerable number of models to find a good fit of the ground state HC5N lines. Again, 
the key parameters are the number density, the temperature and the position of the CCS. 
The position has been fixed between diameters 3.0" and 4.5" after trying several other 
positions in our models. Then, the temperature and density profiles can be established by 
running a grid of models. All the other parameters are then explored to fine tune the shape 
of the line as J up varies. For example, one of the most remarkable differences in the line 
profiles from ground to vibrationally excited HC 5 N is the position of the maximum depth 
of the absorption, its width and the absorption/emission ratio. As shown in Figure 1 and 
table 1, the maximum depth of the absorption for vibrationally excited HC5N happens at 
~ -30 km-s" 1 in the 2 and 3 mm windows (as for the lines in vibrationally excited states of 
HC3N presented in PAR04), whereas for the lines in the vibrational ground state it happens 
at around -40 km-s -1 and is generally wider and less deep. The observational differences 
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have allowed us to determine the kinematic conditions in the CCS. The parameters of the 
best simulation (Fig. 4) are shown on Table 2. The rotational temperature is in the range 
50-70 K. In fact, a simulation with T rot =50 K provides a better fit for lines below J up ~55, 
whereas T rot =70 K provides a better fit for the range J up ~55-80. The simulation shown in 
Fig. 4 is a compromise between the two cases (T rot =60 K). Above J up ~80, the line profiles 
are almost exclusively determined by the SEE region. The error bar for the [HC3N]/[HC5N] 
= 3 ratio in the CCS can again be estimated to be ±1. The same absorption at vlsr <-40 
km-s" 1 reported in section 4.1, and not perfectly modeled here, is seen in some v=0 HC 5 N 
spectra. Absorptions centered at ~ -40 km-s™ 1 have previously been reported in this object 
in molecules such as CO (Herpin et al. 2002), HCN (CER89 and Neri et al. 1992), HCO+ 
(CER89 and Sanchez-Contreras & Sahai, 2004), HC 3 N (CER89, PAR04), and CS (Hajian 
et al. 1995). There is general agreement to assign that absorption to the AGB circumstellar 
envelope. 

Sanchez-Contreras & Sahai (2004) suggest that an elongated structure of ~11" diameter 
with expansion velocity of ~17.5 km-s™ 1 is responsible for part of the emission of the low 
velocity molecular gas in their HCN J=l-0 and HCO + J=l-0 interferometric maps. This 
structure is most probably an extension (due a larger abundance of this molecules 
in cold regions) of the CCS we are modeling here for single antenna data. The difference 
between this larger shell and our CCS is probably due to the different excitation conditions 
of the low-J lines of HCN and HCO + and the high-J lines of HC5N explored in this work. 

On the basis of the extended model, that now includes two gas regions, we can check 
for the presence of longer cyanopolyynes rather easily because only their abundance ratio 
with respect to HC 5 N or HC 3 N is needed to make predictions of their lines and then check 
against the data (assuming the same rotational temperature). This is presented in the 
following subsections. 

4.4. Isotopic HC 5 N and 12 C/ 13 C ratio in the CCS 

We have looked at the possible detection of the 13 C substituted species of HC5N by 
comparing the observations at the frequencies of their rotational lines in the ground vibra- 
tional state with the extended model results considering a 12 C/ 13 C isotopic ratio of 15 in 
the SEE (as derived in section 4.2) and also in the CCS. The simulations indicate that this 
value is an upper limit for the CCS. This is important to later evaluate the HCN and HNC 
abundances from H 13 CN and 13 C (section 4.6). 
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Table 2: Model parameters of the Circumstellar Shell (see figure 3) that provide the best fit 
to the HC 5 N rotational lines in the ground vibrational state (see figure 4). 
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4.5. HC 3 N in the CCS 

The HC3N lines in vibrational states below Iv-j show an emission that exceeds what is 
predicted by the SEE model. Including the CCS component as in section 4.3 provides a first 
approximation to the detected signal (see figure 5) that is in fact quite good for J up below 
~20. However, the HVW component is also visible in the high-J v—0 lines and therefore 
should be included in the model. The high velocity outflow is modeled by just considering a 
gas component with a given excitation temperature (350 K), average column density (2.5-10 17 
cm -2 ), filling factor at a reference frequency (1.3-10 -3 at 90 GHz) and halfwidth velocity 
outflow (150 km-s -1 ). The rather high temperature for this component is confirmed by the 
dominant role that the broad HVW component has at high J's (see CSO lines in figure 5). 
More precise details about the structure of the fast outflow are difficult to obtain with the 
present data. A better constraint on the temperature should be obtained by observing a few 
v—0 HC 3 N lines above 400 GHz (on going). 



4.6. HCN and HNC emission 

At this point of the analysis we have a model that considers the three main regions of 
the gas emission in the millimeter wave spectrum of CRL618: SEE, CCS, and HVW. We 
can try to fit the observed lines from other species by adjusting only their ratio to HC3N 
in the assumption that the rotational temperature is the same. This approximation will be 
quite good for those species arising mainly from the SEE, due to the high density of the 
region (populations are thermalyzed). It could also be useful for the HVW. However, the 
low densities imply that the full analysis for all molecular species will be more complex in 
the CCS. 

In this work we will focus only on HCN and HNC. For opacity reasons, we cannot probe 
the inner SEE and HVW regions using the main isotopologues of both species, specially as 
J increases. Instead we can use H 13 CN and HN 13 C assuming a SEE 12 C/ 13 C value of 15 
(see above) but even for these isotopologues some contribution from the CCS can 
be expected. As a result, our simulations considering only the SEE+HVW (see 
figure 6) for 3 up >2 provide upper limits for the H 13 CN/HC 3 N and HN 13 C/HC 3 N 
ratios in the SEE. The values obtained are about 1/2 and 2/15 respectively which 
translate into HCN/HC3N - 7.5 and HCN/HC3N ~ 2.0. 

To improve the above estimates in the SEE it is better to use vibrationally 
excited states of HCN and HNC since the CCS contribution should be negligible 
for them. The lowest one (t>2=l), at 712 and 477 cm -1 respectively for HCN and HNC, 
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is detected in both molecules (I doublets for all rotational transitions except the forbidden 
J=l-0). We have selected those lines that have been observed and are not badly blended 
with other species to get a better constraint for the above ratios. This analysis (Figure 7) 
results on HCN/HC 3 N ~ 2 and HNC/HC 3 N ~ 1/3 to 1/5. It is thus demonstrated that 
much lower ratios of HCN and HNC with respect to longer cyanopolyynes are found in the 
inner gas of CRL618 (PPNe phase) compared to those found in IRC+10216 (AGB phase). 
See discussion in section 5. 



4.7. Ground state HC 7 N 

Finally, we have studied the HC7N rotational lines within the v=0 state. This molecule 
is clearly detected since all transitions between J up =72 and J up =85 appear above a 3a level 
except for three cases due to a coincidence with a much stronger line (see Table 3 and figure 
8). Note that most of the HC 7 N v—0 signal arises from the CCS. We reach this conclusion 
after running models for which no HC 7 N was present in the CCS, and its ratio to HC 5 N was 
set to be 1/3 in the SEE. In this case, the model predicts some detectable absorption and 
almost no noticeable emission. The fact that HC7N is seen almost exclusively in emission with 
a much larger intensity than in this model, can be easily reproduced if the CCS is included 
with a HC5N/HC7N abundance ratio around 3 and the same rotational temperature for 
both species (see figure 8). The model also predicts that above J up =85 (at an energy level 
of 134 cm -1 ) the signal would be too week to be detected above the 3a level of the survey. 
Similarly, the intensity of lines of HC9N, with a ratio of 1/3 or less with respect to HC7N, is 
estimated to be well below the detection limit of our observations. 



4.8. Vibrationally excited HC 7 N 

The frequencies of the rotational transitions in the lowest energy bending modes of 
HC7N (^13,^14,^15) have become available in a recent work by Bizzochi & Degli Esposti 
(2004). We have searched for those in the lowest vibrational state (0,0,1) but no individual 
detection has been obtained. These lines have been modeled by considering computed ab 
initio values for the (0,0,1), (0,1,0) and (1,0,0) energies: 62, 163 and 280 cur 1 respectively 
(Botschwina et al. 1997). A tentative detection is obtained when performing the addition 
in the vlsr scale of the observed signal at the positions of the 1~ and 1 + 3 up =72 to 79 lines 
(see figure 8). Five lines in the considered J up range were discarded from the summation 
because of blending with stronger lines (see Table 3). 
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Table 3: Detected lines of HC7N toward CRL 618 (laboratory frequencies in GHz are rounded 
to 0.1 MHz). Frequencies in parenthesis correspond to co-added lines in the velocity scale 
for a tentative detection of z^ 15 (see text). 
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5. Comparison with chemical models 



The abundance ratio for consecutive members of the cyanopolyynes family of molecules 
in IRC+ 10216 is ~ 3 (Cernicharo et al. 1987). However, the HCN/HC 3 N ratio in the same 
object is larger than 100. In CRL 618 our results indicate also an abundance ratio near 3 for 
HC3N/HC5N and -6 for HC 5 N/HC 7 N, but for HCN/HC3N it is around 2-3. The difference 
in this ratio at different evolutionary stages of C-rich objects tells us about the chemical 
processing of the gas when the central star evolves toward the white dwarf stage. CER04 
has modeled the role of the strong UV field arising from the hot central star on the chemical 
abundances of the polyacetilenes, cyanopolyynes and carbon clusters. This work shows that 
HCN is quickly photodissociated allowing the production of CN which reacts with C2H2 and 
C2H to produce HC3N and longer cyanopolyynes, thus explaining the low HCN/HC3N ratio 
found in CRL 618. Moreover, this behavior is similar to that found by Cernicharo et al. 
(2001) for C 4 H 2 and C 6 H 2 . The large abundance of HNC with respect to HCN in CRL618 
(HCN/HNC ~ 10 instead of > 200 in IRC+10216, Cernicharo et al. 1987) also points in the 
same direction. It seems that in a short period of time it could be possible to transform an 
important fraction of HCN into HC 3 N and longer cyanide carbon chains. 



6. Summary and future work 

We have extended our previous analysis of the gas shells surrounding the protoplanetary 
nebula CRL618, by studying the pure rotational lines of HC 5 N in its fundamental and the 
lowest 4 vibrationally excited states (first astronomical source in which vibrationally excited 
HC 5 N has been detected), and HC 7 N rotational lines in its fundamental vibrational state. We 
have found that the HC3N/HC5N ratio in the innermost slowly expanding envelope is ~ 3. 
With this ratio, and all the physical parameters previously derived in PAR04, a good match 
to the observed HC5N lines in vibrationally excited states is found. However, the predicted 
lines for the ground vibrational state are far too weak, suggesting that a more extended and 
colder gas shell has to be considered. The physical parameters of this envelope, called cold 
circumstellar shell (CCS) in this work, have been determined. The position of the shell is 
~3.0"-4.5" from the central star and its excitation temperature is only ~60 K (compared to 
250-275 K in the innermost SEE). The column density of the CCS in front of the continuum 
source is negligible compared to that of the SEE. With the CCS derived parameters it is 
found that the observed HC7N lines in the ground vibrational state can be well reproduced 
with a HC5N/HC7N abundance ratio of ~ 3. The isotopic ratio 12 C/ 13 C is found to be ~ 15 
in the SEE from the lines of vibrationally excited HC3N. The same value must be an upper 
limit in the CCS according to the non detection of 13 C substituted isotopologues of HC 5 N. 
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The analysis of the v—0 HC 3 N lines requires to introduce a third component, the 
high velocity outflow, although its structure cannot be precisely known. Finally, the ra- 
tios HCN/HC3N and HNC/HC3N have been derived using the first four rotational transi- 
tions of H 13 CN and HN 13 C, and some lines of vibrationally excited HCN. The result of a 
much lower HCN/HC3N value than in IRC+10216 indicates a processing of HCN into longer 
cyanopolyynes in the protoplanetary nebula stage. 

A first examination of lines from all other molecules detected in the survey shows that 
most of them can be well fitted using the model of the CCS, HVW and SEE gas regions 
described in this paper and in PAR04, so that a good estimate of their abundances in these 
regions can be achieved. Therefore, the next step in our analysis of the CRL618 line survey 
will be to provide a picture of the chemical content in its three main gas components. 

We thank the support of the IRAM-30m staff during the long completion of the line 
survey. CSO operations are supported by U.S. NSF grant AST 22-09008. This work has 
also been supported by Spanish DGES and PNIE grants ESP2002-01627, AYA2002-10113-E 
and AYA2003-02785-E. 

REFERENCES 

Bizzocchi, L., and Degli Esposti, C, 2004, Ap.J. 614, 518. 

Botschwina, P., Horn, M., Markey, K., & Oswald, R., 1997, Mol. Phys., 92, 381 

Bujarrabal, V., Gomez-Gonzalez, J., Bachiller, R., & Martin-Pintado, J., 1988, A&A, 204, 
242. 

Cernicharo, J. Guelin, M., Mart in- Pint ado, J., Pehalver, J., & Mauersberger, M. 1989, A&A, 
222, LI. 

Cernicharo, J., Guelin, M., Menten, K. M., & Walmsley, C. M., 1987, A&A, 181, LI. 

Cernicharo, J., Heras, A. M., Tielens, A.G.G.M., Pardo, J. R., Herpin, F., Guelin, M., 
Waters, L.B.F.M. 2001, Ap.J. 546, L123. 

Cernicharo, J., 2004, Ap.J., 608, L41. 

Cernicharo, J., & Pardo, J. R., 2005, A&A, In preparation. 
Hajian, A.R., Phillips, J. A., Terzian, Y., 1995, ApJ, 446, 244. 



- 15 - 

Herpin, F., Goicoechea, J.R., Pardo, J.R., Cernicharo, J. 2002, Ap.J., 577, 961 

Kahane, C, Cernicharo, J., Gomez- Gonzalez, J., Guelin, M., 1992, A&A, 256/1, 235. 

Kroto, H.W., Kirby, C, Walton, D.R.M., Avery, L.W., Broten, N.W., MacLeod, J.M., Oka, 
T., 1978, ApJ, 219, L133. 

Little, L.T., MacDonald, G.H., Riley, P.W, Matheson, D.N., 1978, MNRAS, 183, 45R 

Neri, R., Garcia-Burillo, S., Guelin, M., Cernicharo, J., Guilloteau, S., Lucas, R., 1992, 
A&A, 262, 544. 

Pardo, J.R., Cernicharo, J., Goicoechea, J.R., Phillips, T.G., 2004, Ap.J., 615, 495. 
Sanchez-Contreras, C, & Sahai, R., 2004, Ap.J., 602, 960. 

Snell, R.L., Schloerb, F.P, Young, J.S., Hjalmarson, A., Friberg, P., 1981, ApJ, 244, 45S. 

Yamada, K.M.T., Degli Esposti, C, Botschwina, P., Frster, P., Bizzocchi, L., Thorwirth, S., 
Winnewisser, G., 2004, A&A, 425, 767. 



This preprint was prepared with the AAS IATgX macros v5.2. 



-16- 



J up > 



35 



57 



84 



vib. 
state 



Q ~ 93.5 GHz) (v ~ 152 GHz) {v ~ 224 GHz) ^ 




i/ n le 



^io le 



u u If 



l/ 10 If 



60 -40 -20 -60 -40 -20 -60 -40 -20 




2v 1] 



-100 100 -100 100 -100 100 
V LSR ( km O 

J im — > 31 35 39 

up state 
(v - 83.0 GHz) (v - 93.8 GHz) (v ~ 104.5 GHz) 




V LSR ( km S V ) 



Fig. 1. — Sample of detected lines in vibrationally excited states of HC5N toward CRL 
618. Model results, shifted by 0.2 F C (F c : Continuum Flux) for all vibrational 
states except for 3u u (0.1 F c shift), are shown at the velocity resolution of the data 
(Av = c - ^, with Au ~1 MHz for all spectra) for a HC3N/HC5N value of 3. All other 
physical conditions are as determined in PAR04 in the SEE (the only region considered). The 
extra feature that appears in the 2uu J up =35 spectrum is a blending of the {351 recombination 
line of atomic hydrogen plus the J=83-82 line of v=0 HC7N. 
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Fig. 2. — Model predictions considering a 12 C/ 13 C isotopic ratio of 15 in the slowly expanding 
envelope of CRL 618 (shifted by 0.2F C ) compared to the observations for a sample of lines 
belonging to single 13 C substituted HC3N in the v 7 state. The velocity resolutions of 
model and data are the same (Av — c ■ with Av ~1 MHz for all spectra and 
v ~J up -9.1 GHz). Two weak features, between the HCC 13 CN and HC 13 CCN vl l~ lines, 
clearly seen in the spectra with the best signal-to- noise ratio, correspond to the uQ 1 + lines 
of both species. They have also been included in the model. All other extra features that 
appear in some spectra belong to vibrationally excited states of HC 3 N. 
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Fig. 3. — Model of CRL 618 including a cold circumstellar shell first introduced in our work 
to reproduce the observed rotational lines of HC5N in its ground vibrational state. Details 
in section 4.3. 
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Fig. 4. — (Left) Observed rotational lines of HC 5 N in the ground vibrational state compared 
with the SEE model assuming a HC3N/HC5N value of 3 (dotted histogram, shifted by 
1.0-F C ) and with the expanded model (solid histogram, shifted by 1.0 F C ), that 
includes the CCS, presented in this paper (Sect. 4, Fig. 3 and Table 2). Both 
models are basically coincident for 3 up above ~ 80. (Right) Results of the CRL 618 
model including the SEE and the CCS regions (solid histogram shifted by 0.2-F c ) for 
a sample of lines of 13 C substituted HC5N, compared to the observations. The simulations 
consider a 12 C/ 13 C value of 15 in both regions. Velocity resolution of data and model 
are the same for each plot, as explained in previous figures. 
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Fig. 5. — Observed rotational lines of HC3N in the ground vibrational state compared 
with the SEE+CCS model assuming a HC3N/HC5N value of 3, plus a HVW component 
according to section 4.5 (solid histogram shifted by 1.0 F C ). The model for the SEE 
region alone is also plotted for comparison (dotted line, shifted by -1.0 F C ). Some 
other lines appear in the data and are also modeled. Most of them are rotational transitions 
within different vibrational states of HC3N, HC5N or their isotopologues, but there are some 
belonging to other species the discussion of which we leave for a separate work: H 2 CO and 
C5H in J up =16, and C-C3H2 in J np =35. All observations have been carried out with the 
IRAM-30m telescope except J up =31,35 (CSO). Velocity resolutions of data and model 
are the same as explained in previous figures. 
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Fig. 6.— Observed rotational lines of H 13 CN and HN 13 C compared with the SEE+HVW 
model (shifted by 2.0 F C ) assuming H 13 CN/HC 3 N and HN 13 C/HC 3 N values of 1/2 and 
2/15 respectively. All observations have been carried out with the IRAM-30m telescope 
except J up =4 (CSO). The very low excitation temperatures of J up >=1,2 levels, compared 
to those in the SEE and HVW regions, results in the the poor agreement between data and 
model for the lines arising from those levels. Nevertheless, the other lines provide a useful 
upper limit (the CCS contribution should be important) about the HCN and HNC 
abundances in the SEE and HVW regions. Velocity resolutions of model and data are 
the same for the different lines, as explained in previous figures. 
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Fig. 7. — Observed rotational lines of HCN and HNC in the v=2 vibrational state 
compared with the SEE model (shifted by 2.0 F C ) assuming HCN/HC 3 N and HNC/HC 3 N 
values of 2 and 1/5 respectively. All observations have been carried out with the IRAM-30m 
telescope except J up =4 (CSO). Velocity resolutions of data and model are the same 
as in previous figures. 
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Fig. 8. — HC 7 N (v—0) lines observed toward CRL618 compared with model results including 
the SEE and CCS regions (solid histogram shifted by 0.075 F c ) and the CCS region alone 
(dotted line also shifted by 0.075 F C ). The considered value of HC 5 N/HC 7 N is 3, 

which provides a very good fit. The spectrum in the lower-right corner box is the addition 
of the z/ 15 HC7N rotational lines listed in table 3 compared to the model results for both the 
CCS and the SEE. A 4.0 scaling factor has been introduced to keep the same scale of all 
other spectra. 



